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' Abstract. We present if-band number counts for the faint galaxies in the Calar Alto Deep Imaging Survey 

(CADIS). We covered 4 CADIS fields, a total area of 0.2 deg 2 , in the broad band filters B, R and K. We detect 
about 4000 galaxies in the K-band images, with a completeness limit of K = 19.75 mag, and derive the Tf-band 
galaxy number counts in the range of 14.25 < K < 19.75 mag. This is the largest medium deep if-band survey 
to date in this magnitude range. The B- and 7?-band number counts are also derived, down to completeness 
limits of B — 24.75 mag and R = 23.25 mag. The if-selected galaxies in this magnitude range are of particular 
interest, since some medium deep near-infrared surveys have identified breaks of both the slope of the if-band 
number counts and the mean B — K color at K = 17 ~ 18 mag. There is, however, a significant disagreement 
in the if-band number counts among the existing surveys. Our large near-infrared selected galaxy sample allows 
us to establish the presence of a clear break in the slope at K — 17.0 mag from dlogiV/dm = 0.64 at brighter 
magnitudes to dlogiV/dm = 0.36 at the fainter end. We construct no-evolution and passive evolution models, 
and find that the passive evolution model can simultaneously fit the B-, R- and Jf-band number counts well. The 
B — K colors show a clear trend to bluer colors for K > 18 mag. We also find that most of the K = 18 — 20 mag 
galaxies have a B — K color bluer than the prediction of a no-evolution model for an L» Sbc galaxy, implying 
either significant evolution, even for massive galaxies, or the existence of an extra population of small galaxies. 
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1. Introduction 

Many galaxy surveys have shown that galaxies were un- 
dergoing a significant evolution at intermediate redshifts 
(0. 5 < z < 1.0). Cowie et al. (|1996|) and Ellis fll997| , et 
al. 1996] ) suggested that an extra galaxy population in this 
redshift range is responsible for the excess of the galaxy 
number counts at faint magnitudes. With additional mor- 
phological classification for the faint galaxies, this extra 
population was identified as irregular and peculiar galax- 
ies with strong star formation at intermediate redshifts 
(Glazebrook et a l. |1995| , Abr aham et al. |1996| , van den 



Bergh et al. 1996, Huang et al. 1998). The massive ellipti- 
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cal and giant spiral galaxies, however, show no significant 
evidence of evolution since z = 1.0 (Lilly et al. 1995 ). 



As an indication that the faint blue galaxy excess has 
a significant contribution from galaxies at high redshift 
Cowie et al. (1997) find that at least some of their faint B- 
selected galaxies are located in a high- redshift (1 < z < 2) 
tail on the redshift distribution. Dickinson ( 2000] ) con- 
ducted a near-infrared study on HDF-N galaxies at using 
the NICMOS camera on HST and found far fewer high 
luminosity galaxies in the 1.4 < z < 2 range. However, 
the NICMOS data cover only a small field, and may not 



be a representative sample. Steidel et al. (1996d999) have 
identified a large number of galaxies at higher redshifts 
(z ~ 3 — 4) using the Lyman-break technique. However, 
these galaxies do not contribute significantly to the faint 
blue galaxy number counts at B < 26 mag. For z > 3, the 
B band covers wavelengths below the redshifted Lyman- 
alpha (A rost = 1216 A) line, where the rest-frame UV light 
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Table 1. Field coordinates 



Field 


a(2000) 


5(2000) 


area [arcmin 2 ] 


01H 


01:47:33.3 


02:19:55 


175 


09H 


09:13:47.5 


46:14:20 


127 


16H 


16:24:32.3 


55:44:32 


186 


23H 


23:15:46.9 


11:17:00 


188 



is strongly affected by absorption in the Ly-a forest as 
well as dust. 

Comparing the galaxy luminosity function at differ- 
ent redshifts is the most direct way of understanding the 
galaxy evolution. However existing samples are not large 
enough and the galaxy luminosity function at interme- 
diate redshifts is not well determined (Lilly et al. 1995, 
Cowie et al. 1996, Ellis et al. 1996), especially at lower lu- 
minosities. A wide field medium deep survey is required to 
obtain larger numbers of low luminosity galaxies, to better 
constrain the luminosity function at intermediate redshift. 
While large-format CCDs have been in use for some time 
now, sufficiently large near-infrared arrays have only re- 
cently become available. 

There are three significant advantages to studying 
galaxy evolution in the near-infrared K band instead of in 
the optical bands (Cowie et al. 1994, Huang et al. 1997). 
First, the cosmological k-correction in the K-hanA is much 
smaller and better determined than the k-corrections in 
the optical bands. Second, the shorter wavelength optical 
bands sample restframe blue or UV light, and are thus 
very sensitive to the star-formation history in galaxies. 
Third, the if -band light is a close tracer of the stellar 
mass in galaxies, and less sensitive to either the star for- 
mation history or the galaxy spectral type. A if-selected 
sample of galaxies is thus more representative of the true 
distribution of galaxies. 

A number of medium deep and deep if-band galaxy 
surveys have been conducted by various groups (Gardner 
et al. |1993[ S ongai l a et a l. |1994| , Djorgo vski e t al. |l995 , 
Sarac co et al. |l997j |l999j M ineza ki et al. |1998| , Bershady 
et al. 1998, Moustakas et al. 1998 and Kiimmel & Wagner 



200C). There is some disagreement between the if-band 
number counts of these surveys at faint magnitudes. Still, 
these surveys do show evidence for two critical phenomena 
which are important in understanding galaxy evolution: 
the mean B — K color becomes bluer at faint magnitudes, 
and the slope of the K-b&nd number counts shows a break 
in the 17 < K < 19 magnitude range. The Hawaii if-band 
galaxy survey (Songaila et al. 1994) shows that galaxies 
in the range of 16 < K < 20 magnitude are located at 
0.3 < z < 1.3, which spans an important epoch for galaxy 
evolution. 

To better determine the if-band number counts in 
this magnitude range, a if-band survey with much larger 
coverage is essential. Here, we use the Calar A lto D eep 
Imaging Survey (CADIS, Meisenheimer et al. 1998 ) to 
derive improved if-band number counts in the range 
14.25 < K < 19.75 mag, as well as the corresponding 



optical number counts in the B and R bands. CADIS 
was originally tailored to search for Ly-cv emitters at high 
redshift with Fabry-Perot Interferometer (FPI) imaging 
(Thommes et al. 1998). However, the CADIS fields are also 
observed in three broad band filters (B, R, J or K) and 
up to 13 medium band filters with A/AA = 25 — 50. This 
multi-color database is also being used to study galaxy 
populations and their evolution at intermediate redshifts 



through their photometric redshifts (Wolf 1998), as well as 
the detection of unusual objects like active galactic nuclei 



(Wolf et al . 1999) and e xtrem ely red galaxies (Thompson 
et al. |1999|) . Fried et al. (|200l|) present initial results from 
CADIS on the luminosity function of optically selected 
galaxies and their evolution. 

In this paper, we present the optical and near-infrared 
galaxy number counts, as well as the colors of the K- 
selected galaxies. The four CADIS fields used for this 
study cover a total of 0.2 deg 2 down to a completeness 
limit of K = 19.75 mag. Our coverage is four tim es th e 
area of the ESO K-b&nd survey (Saracco et al. 1997 ), 



which is the largest medium deep i^-band survey to date. 
The paper is organized as follows: The observation, data 
reduction and photometry are described in Section 2; 
the B-, R-, and if-band number counts are presented in 
Section 3, along with the colors of the if-selected galaxies. 
In Section 4 we present a galaxy model which reproduces 
the observed number counts, and Section 5 summarizes 
our results. 



2. Observation and photometry 

2.1. Observation and data reduction 

A complete description of the observation and data reduc- 
tion for CADIS will be given elsewhere (Meisenheimer et 
al. in prep.), here we give only a brief description. The 
CADIS fields are selected at high galactic latitudes on 
the northern and equatorial sky. The two primary CADIS 
field selection criteria are: (1) no bright star (R < 16 mag) 
within a field of 12'xl2', and (2) a local minimum on the 
IRAS 100 [im maps, with an absolute surface brightness 
of < 2MJy/sr. The field positions and areas covered in 
this project are given in Table 0. 

The optical observations were obtained with focal re- 
ducing CCD cameras on the 2.2m and the 3.5m tele- 
scope on Calar Alto Spain. The Calar Alto Faint Object 
Spectrograph (CAFOS) covers a 200 arcmin 2 field at 0.53 
arcsec per pixel on the 2.2 m telescope, while the Multi- 
Object Spectrograph for Calar Alto (MOSCA) covers a 
150 arcmin 2 at 0.33 arcsec per pixel on the 3.5 m telescope. 



1 The CADIS I filter used here is a narrow band filter at 
A = 815 nm (<5A = 25 nm), and used only for the star-galaxy 
classification. The object detection program was not run on 
the I-band images. 

2 The Calar Alto 2.2 m and 3.5 m telescopes. 

3 The integration time in seconds. 

4 Mgo is the magnitude where 90% of the galaxies are de- 
tected in our Monte Carlo test. 
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B 






R 






I 1 




K 




Field 


Tel 2 


Int 3 


M 90 | 


Tel 


Int 




Tel 


Int 


Tel 


Int 


M 90 


01H 


2.2 


5050 


23.7 


2.2 


2800 


22.8 


2.2 


4000 


3.5 


9648 


19.4 


09H 


3.5 


4000 


24.9 


2.2 


2800 


23.4 


3.5 


7000 


3.5 


11280 


18.8 


16H 


2.2 


6200 


24.6 


2.2 


2500 


23.7 


2.2 


13620 


3.5 


9067 


19.3 


23H 


2.2 


6230 


23.4 


2.2 


2500 


23.0 


2.2 


5000 


3.5 


6475 


19.6 



Only images with seeing better than 1'.'8 are included in 
the CADIS database. 

The if-band data were obtained with the Omega- 
Prime near-infrared camera (Bizenberger et al. 1998] ) on 
the 3.5 m telescope at Calar Alto. Omega-Prime has a 
1024x1024 pixel HgCdTe array which, at 0.396 arcsec 
per pixel, covers a 6.8' x 6.8' field of view. This camera 
has no cold pupil stop, which would effectively block ther- 
mal emission from the warm surface around the pr imary 
mirror. We used a K'-Sitei (Wainscoat & Cowie 1992 ) 
instead of the standard -fT-band to reduce the thermal 
background. In order to construct a larger field of view, 
matching that of the optical cameras, we observed the 
CADIS-fields in a 2 x 2 mosaic. After trimming the high 
noise edges of the stacked mosaics, the final coverage for 
the if-band images varies from 127 arcmin 2 to 188 arcmin 2 
(see Tables |l| and f|. The K band image quality criterion 
is a seeing less than 1'.'5. 

The optical images were initially processed using the 
Munich Image Data Analysis System (MIDAS) software 
package, while the near-infrared images were processed 
using the Interactive Data Language (IDL) or IRAF^j. 
Thompson et al. (1999) gives a discussion of the steps of 
the standard reduction to get the final, stacked images in 
the optical and the near-infrared. A detailed discussion 
on the dctectability and photometry is given in the next 
section. 

Photometric calibrations for all of the CADIS data are 
determined on the Vega scale. Each CADIS field has one 
or two relatively bright stars for which we have established 
accurate spectro-photometric calibrations between 350 nm 
and 950 nm. The magnitudes of these tertiary standards 
are then obtained by integrating their spectral energy dis- 



tributions within the CADIS filter profiles. Wolf ( |1998[ ) 
gives a detailed description of the CADIS magnitude sys- 
tem and its calibration. For the A"-band data, we use the 



UKIRT faint standard stars (Casali & Hawarden 1992) 
for flux calibration. The total exposure times and detec- 
tion limits for the optical and infrared data are given in 
Table [§ 

We use the CADIS broad band B, R, and K data for 
the number counts, and add the medium band / (<5A = 
25 nm, centered at 815 nm) data in order to facilitate star- 



J IRAF is distributed by the National Optical Astronomy 
Observatories, which are operated by the Association of 
Universities for Research in Astronomy, Inc., under cooper- 
ative agreement with the National Science Foundation 



galaxy separation. The CADIS B and R filters are very 
close to the Johnson- Cousins system (Bj and Rkc in the 
equations below). The transformations between the two 
systems were derived using synthet ic ph otometry of the 
galaxy templates of Kinney et al. ( 1996 ) in the redshift 
range < z < 1.5. The transformation equations are: 

Bj - Bcadis = 0.082 + 0.121(5 - i?) CAD is (1) 

Rkc - i? C ADis = -0.027 - 0.034(5 - i?) C ADis (2) 

2.2. Object detection and photometry 

Object catalogs were created using the S Extractor image 
analysis package (Bertin & Arnouts 1996 ). Before running 
the detection program, each image was convolved with a 
Gaussian PSF with a width equal to the full width at half 
maximum (FWHM) of the seeing. For detection, we re- 
quired at least 5 connected pixels, each with a signal larger 
than three times the background noise. Faint galaxy de- 
tectability is determined using Monte Carlo simulations. 
Several bright galaxies with different morphological types 
were rcscalcd then added back into the original image be- 
fore running the detection program again. The result is 
the percentage of galaxies recovered as a function of ap- 
parent magnitude. The magnitudes at which 90% of the 
test galaxies are recovered are listed in Table |[ 

We used the photometry package MPIAPHOT 



(Meisenheimer & Roser 1986), which was developed to 
measure accurate spectral energy distributions (SED) 
from images with variations in seeing. In MPIAPHOT, 
the pixel values within a circular aperture are evaluated 
using a Gaussian weighting function. The e-folding width 
s w of the Gaussian weighting function is set individually 
for each frame i to 



'off 



(3) 



Here s, ; denotes the e-folding width of the seeing where x = 
a,b corresponds to the major/minor axis of the empirical 
seeing profile and s c g denotes the e-folding width of the 
effective PSF. The effective seeing is chosen as a constant 
for the dataset of every CADIS field, and compensates 
for variations in the seeing between images obtained at 
different times and through different filters. Since most of 
our images have a seeing of about 1'.'5 FWHM (s, = 0'.'90) 
and the signal to noise ratio of the photometry has its 
maximum at s e g ~ V%Si, the effective seeing is set to 
s c fi = 1'.'38 for all data presented here. 
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Fig. 1. The color-color diagram for the if-selected sample. 
In the diagram, the dots represent the point-like objects 
and the crosses represent the extended objects. The solid 
line represents the empir ical c olor criterion f or sta r-galaxy 
separation (Cowie et al. 1994, Huang et al. 1997). 



The weighted pixel values are summed within an aper- 
ture of 4" diameter, and then corrected to a total mag- 
nitude. By default MPIAPHOT applies corrections calcu- 
lated by integrating the flux of point-like objects to infi- 
nite radius. The corrected magnitude thus equals the total 
magnitude for point sources, but underestimates the total 
magnitude of extended galaxies. To correct for this effect, 
we apply another correction for galaxies. The corrections 
are determined from simulations containing galaxies which 
match the range of apparent sizes seen in the survey data. 
We find that the corrections for resolved galaxies are al- 
most independent of the bulge to disk ratios and inclina- 
tions. The correction can be well described as a function of 
the measured size of the galaxies (i.e. the second moment 
along its major and minor axis) out to about 3" FWHM. 
The derived correction formula has been applied to the 
science data, in order to convert CADIS magnitudes to 
total magnitudes. 

2.3. Star-Galaxy Separation 

Conventional star-galaxy separation is based on compar- 
ing the morphology of objects with the morphology of 
point-like objects. Since the if-images have the largest 
field size, there are regions where only information in K is 
available, and there we used this conventional method for 
star-galaxy separation. However, this method fails to iden- 
tify compact galaxies with unresolved morphology. While 
those objects are rare at bright magnitudes, this popu- 
lation amounts to ~ 20% of all galaxies at R = 23 mag. 
To identify this population we use colors to separate stars 



and galaxies. In three fields presented here (Olh, 09h and 
16h field) the full CADIS dataset in broad bands and up 
to 13 medium band filters is available. In those fields the 
broad wavelength coverage from 0.4 /im to 2.2 /im allows 
a classification into stars, galaxies and QSOs based on the 
object's SED. This is done by comparing the measured 
color indices of every object with the expected color in- 
dices calculated from template spectra of the different sub- 



classes stars, galaxies and QSOs (Wolf 1998). For galax- 
ies and QSOs the theoretical color indices are calculated 
with a broad range of redshifts (0 < z < 2 for galax- 
ies and < z < 6 for QSOs), the result of the classi- 
fication includes a redshift estimate with an accura cy of 
typically Sz = 0.03 for galaxies (Fried et al. |2001 ). For 
galaxie s, we use the set of templates published by Kinney 
et al. ( 1996 ) whose SEDs range from an elliptical to an 
extreme starburst galaxy. This further allows a classifica- 
tion of the spectral type of galaxies found in the survey. 
For the 23h-field, for which we do not yet have the full 
color information, color-based classification was done us- 
ing the color-color dia gram I — K vs. B — I (Cowie et 
al. 1994, Huang et al. [1997 ). To demonstrate this tech- 
nique we plot in Fig. |l| the color-color diagram I — K vs. 
B — I of pointlike (dots) and extended (crosses) objects 
for the -fT-selected sample. In the if-selected sample stars 
and galaxies are well separated by the empirical relation 



B - I = 2.5 * (I - K) - 2 



(4) 



with only a few galaxies crossing the line. 

To summarize the criterion for the final star-galaxy 
separation: an object is classified as a galaxy if either it 
is morphologically identified as a galaxy, or it is identified 
as a galaxy on the basis of its colors. 

3. Optical and near-infrared galaxy number counts 

The galaxy number counts in the three filters are listed in 
Tabl es |3f ]5[. The number counts in B and R are plotted in 
Figs. |2| and p[ respectively. We also compare our number 
counts to those of other large scale surveys with areas on 
a one square degree scale, like Arnout s et a l. ( 1997 ), EIS 
(ESO Imagi ng Su rvey, Prandoni et al. 1999 ) and Kiimmel 
& Wagne r ( 2001 ) and th e deep surveys from Metcalfe et 
al. fll995| ), Smailet al. fll995| ), the Hub ble Deep Field 
(Williams et al. |1996[ ), Hog g et al. ( |1997[ ) and the NTT 
Deep Field (Arnouts et al. 1997) In both filters there is 
a very good agreement between the CADIS data and the 
average literature counts over the entire magnitude range. 
Apart from serving as input for the models presented in 
the next section the agreement of our optical data with 
the literature confirm that the CADIS-fields contain the 
typical mix of galaxies, as it is expected for the wide field 
of view and the four lines of sight. 

The CADIS if-band counts are plotted in Fig. |. We 
also plot for comparison t he re sults from other if -band 
surveys ( Gard ner e t al. 1993 Djo rgovsk i et a l. 1995 , 
Gard ner et al. |l996| , Sarac co et al. |l997j , |l999| , Huang 
et al. h997l Minezaki et al. Il998|, Moustakas et al. 11998, 
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Fig. 2. The B-band galaxy number counts. The filled cir- 
cles are the number counts obtained in CADIS. The other 
symbols are number counts from Metcalfe et al. (1995) 



(Met), Williams et al. ( [1996 ) (Wi ll), Arnouts et al. ( p.997| ) 
(Arn97), Prandoni et al. (1999) (Pra) and Arnouts et 
al. ( 1999| ) (Arn99). The upper pair of lines show the pas- 
sive evolution models, while the lower pair of lines show 
the no evolution models. In both cases, the solid lines are 
for go = 0.05 while the dashed lines are for go — 0.5 mod- 
els. 



Bersh ady et al. |l99jj, Szokoly et al., (1999), McCracken 
et al. 2000 and Kummel & Wagner 2000). Because of 
the wide range in magnitudes covered by these surveys 
(14 < K < 24 mag), and the relatively small differences 
between them, we plot the if-band number counts in three 
magnitude ranges for better resolution. 
At magnitudes brighter than K = 17.5 mag, substantial 
-fT-band surveys with areas on the sq uare-degree sca le wer e 
carried out by Gardner et al. ( 1996 ), Huang et al. ( 1997| ), 
Szokoly et al. (1999) and Kummel & Wagner (2000). 
Therefore, the number counts in this magnitude range are 
well defined. As shown in the up-left panel of Fig. |J, the 
CADIS if-band number counts in this magnitude range 
are in good agreement with previous number counts, but 
seem to support only the highest counts in the range 
15.5 < K < 17.5 mag 

In the magnitude range of 17.5 < K < 19.5 mag (up- 
per right panel of Fig. ||) only the surveys of Minezaki 
et al. fll998[ 181arcmin 2 ) and Saracco et al. ( |1997 , 
170 arcmin ) have a coverage comparable to CADIS. With 
an area of 675 arcmin 2 CADIS is the largest survey in this 
magnitude range. The CADIS counts are slightly lower 
than those of Minezaki et al. (1998) at the faint end, 
but 50% higher than the counts of Saracco et al. (1997) 
across the whole brightness interval. Because of the small 
statistical uncertainties due to the large 0.2 deg 2 area 
of the four CADIS fields, our data provide a more ro- 
bust measurement of the number counts in the range 
16.5 < K < 19.5 mag. The largest remaining discrepancies 
in the number counts are for K > 19.5 mag, i.e. beyond 
the limit of the CADIS counts. At these faint magnitudes 
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Fig. 3. The i?-band galaxy number counts. The filled cir- 
cles are the number counts obtained in CADIS. The other 
symbols are number counts from Metcalfe et al. (1995) 



(Met), Smail et al. fll995|) (Sm ail), Arnouts et al. dl997| ) 
(Arn97), Hog g et al. fll997| ) (Hogg) and Kummel & 
Wagner ( |200 1[ ) (Kue) . The models for the i?-band number 
counts are plotted as in Fig. |^. 



the scatter in the counts reach up to a factor of three (at 
K ~ 21.0, see lower left panel in Fig. |). From the CADIS 
counts above the onset of this scatter at K > 18 mag, we 
conclude that CADIS favours the largest values obtained 
for the number counts at K > 19 mag, consistent with 
those obtained by Moustakas et al. (1997). 

Magnitude errors, field-to-field variation and the 
Poisson noise are prime contributors to the uncertain- 
ties in the number counts. Fig. |4| shows in the lower right 
panel the number counts for the individual CADIS-fields. 
The differential counts for the individual fields give consis- 
tent results and therefore we exclude the possibility that 
the average number counts of all fields are affected by 
field to field variations. At K > 18 mag those variations 
can be significant only when the sky coverage of a survey 
is <C 100 arcmin 2 . Accordingly, the field-to-field variation 
and the Poisson noise may have significant contribution 
to the scatter of the if-band number counts among the 
deep surveys which cover only a small area on the order 
of 1 arcmin 2 . 

In most cases, the uncertainties in the number counts 
result from systematic errors of the magnitude scale. As 
it is outlined in Huang et al. (1997), a random magni- 
tude error can cause the slope of the number counts to 
become steeper, while a systematic magnitude error can 
cause either a shift of number counts in the magnitude 
direction or change of the slope. The absolute CADIS- 
photometry (with respect to Vega) is accurate to 0.1 mag 
(Meisenheimer et al. in prep.). Taking into account those 
photometric errors leads to perfect agreement between our 
counts and the Minezaki et al. (1998) data, but even then 
the CADIS counts are still incompatible with the results 
of Saracco et al. (1997). 
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Fig. 4. The if -band galaxy number counts. The filled circles are the number counts obtained in CADIS. The other 
symbols are number co unts f rom Gardner et al. (1993) (Gar93), Djorgovski et al. 1995 ( Djo), Gardner et al. (1996) 



(G ar96) Huang et al. (|l997|) (H uang9 7), Saracco et al. (|l997|) (Sa r97), Minezaki et al. |l998|) (Min), Moustakas et 
al. |l998| (M ou), Bershady et al. (|l998|) (Ber), S zokoly et al., (|l999|) (Szo), Saracco et al. ( |l999|) (Sar99), McCracken 
et al. ( 2000|) (McCra) and Kiimmel & Wagner ( 2000| ) (Kue). The models for the if -band number counts are plotted 
as in Fig. ti. 



Gardner et al. ( 1993| ) reported finding a break in the Table 3. £?-band number counts 



logarithmic slope of the if -band counts at K ~ 17, where 
we have good statistics. To test for this, we fit a power law 
with the form N(m) — a* iQ b * m to our counts in various 
magnitude ranges. We find a clear break in d\ogN/dm 
from a bright-end slope of b = 0.64 to a faint-end slope of 
0.36, with the break in the interval 16.5 < K < 17.0 mag. 
The slope at the bright end is comparable to that found in 



other surveys (see Kiimmel & Wagner 2000). The slope at 



the faint end is somewhat steeper than values reported in 



the smaller surveys of Gardner et al. (1993) and Minezaki 



et al. (1998), who give slopes around 0.24. 



4. Modeling the Number Counts and Galaxy 
Population 

Models of galaxy number counts are primarily constrained 
by assumptions on the cosmological geometry, galaxy 



B mag 


Log(N)0 


Log(N-cr) 


Log(N+a) 


17.0 


1.07 


0.54 


1.30 


18.0 


1.07 


0.54 


1.30 


18.5 


1.37 


1.07 


1.55 


19.0 


1.62 


1.41 


1.76 


19.5 


2.03 


1.91 


2.12 


20.0 


1.97 


1.85 


2.07 


20.5 


2.44 


2.37 


2.50 


21.0 


2.51 


2.45 


2.57 


21.5 


2.83 


2.79 


2.87 


22.0 


3.07 


3.04 


3.10 


22.5 


3.35 


3.32 


3.37 


23.0 


3.60 


3.58 


3.61 


23.5 


3.85 


3.83 


3.86 


24.0 


4.03 


4.02 


4.05 


24.5 


4.25 


4.24 


4.27 



The number counts are in units of N/0.5 mag/deg 2 
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Table 4. i?-band number counts 



23.25 < B < 24.75 



R mag 


Log(N) 


Log(N-cr) 


Log(N+<r) 


16.0 


0.69 


0.15 


0.99 


16.5 


1.39 


1.13 


1.55 


17.0 


1.17 


0.79 


1.36 


17.5 


1.39 


1.13 


1.55 


18.0 


1.87 


1.74 


1.96 


18.5 


2.10 


2.01 


2.18 


19.0 


2.44 


2.38 


2.50 


19.5 


2.68 


2.63 


2.72 


20.0 


2.86 


2.82 


2.89 


20.5 


3.11 


3.08 


3.13 


21.0 


3.27 


3.25 


3.30 


21.5 


3.47 


3.45 


3.49 


22.0 


3.60 


3.59 


3.62 


22.5 


3.78 


3.77 


3.80 


23.0 


4.05 


4.04 


4.06 



Table 5. K-hand number counts 



K mag Log(N) Log(N-a) Log(N+a) 



14.0 


1.21 


0.83 


1.40 


14.5 


1.68 


1.51 


1.81 


15.0 


2.01 


1.90 


2.10 


15.5 


2.32 


2.25 


2.39 


16.0 


2.59 


2.51 


2.62 


16.5 


2.94 


2.91 


2.97 


17.0 


3.12 


3.09 


3.14 


17.5 


3.34 


3.32 


3.36 


18.0 


3.50 


3.48 


3.51 


18.5 


3.69 


3.68 


3.70 


19.0 


3.86 


3.85 


3.87 


19.5 


4.05 


4.04 


4.06 



spectral evolut ion, and the galaxy luminosity function (see 
Gardner 1995). Theoretical uncertainties on the number 



count models come mainly from modeling the spectral evo- 
lution of galaxies, especially in the infrared bands. Many 
galaxy spectral synthesis models currently exist, ranging 
from simple stellar population models assuming an in- 
stantaneous burst of star formation to models including 
gas, recycling of metals, and nebular emission (Bruzual & 
Chariot |l993j , Bert elli et al. |1994| , Worthey |1994| , Fioc & 
Rocca-Volmerange 1997). 

Most models make predictions consistent with each 



other at visible wavelengths (Chariot 199£). There is, how- 
ever, a significant scatter between model predictions at 
near-infrared wavelengths. This is due primarily to poor 
knowledge of the atmospheric parameters for cool stars 
and a substantial contribution from the nebular emission 



(Fioc & Rocca-Volmerange 1997). Many number count 
models (Metcalfe et al. |1996j McCracken et al. |2000D 
use simple stellar population models convolved wi th var - 
ious star formation histories (Bruzual & Chariot 1993|) , 
which over-predict the number counts at bright magni- 
tudes at both visible and near-infrared wavelengths. We 
use PEGASE (Projet d'Etude des Galaxies par Synthese 



30000 - 



20000 - 



10000 



J 



E/S0 



Spiral Sd/Irr 



21.75 < R < 23.25 



30000 - 
20000 - 
10000 - 

— 



E/S0 



Spiral Sd/Irr 



18.25 < K < 19.75 



30000 



20000 



10000 



E/SO Spiral Sd/Irr 

Fig. 5. The observed spectral type mix of galaxies 
(shaded) in the magnitude ranges close to the respective 
completeness limit in comparison to the predicted type- 
mix of the no-evolution and the passive-evolution model 
(the white and black bins, respectively). 



Evolutive), a new spectrophotometric evolution model 
which includes gas and dust components, nebular emis- 
sion, and a better description of galaxy SEDs in the near- 
infrared (Fioc & Rocca-Volmerange 1997). Our number 
count model constructed using PEGASE is able to fit the 
observed number counts at bright magnitudes well (Fioc 
& Rocca-Volmerange 1997). 

We adopt the passive evolution and no-evolution num- 
ber counts models from Huang (1997). The cosmological k- 
correction and evolution e-correction are calculated using 
the PEGASE model. We plot model predictions against 
the observed number counts in Figs. [l]-[|. Huang (1997) 
has shown that these models correctly predict the sur- 
face density and color distribution of local galaxy popu- 
lations, consistent with Fioc & Rocca-Volmerange (1997). 
The model Hubble types refer to the spectral type and 
not the morphological type, although they are correlated. 
The spectral types are defined by the e-folding time of the 
star formation rate. Most models use an exponentially de- 
creasing star formation rate for all spectral types except 
Irr galaxies, which are so blue that the star formation rate 
must have been constant or increasing over the ~4.5Gyr 
life of the galaxy. 

We can see in Figs. that evolution has the strongest 
effect on the models, while the geometry (qo = 0.05 and 
0.5) does not make a large contribution at brighter mag- 
nitudes (B < 24 mag, R < 23 mag and K < 20 mag). We 
adopt a high normalization of (/)„ = 0.02 h 3 Mpc~ 3 . Only 
the passive evolution model is able to simultaneously fit 
the B-, R- and X-band number counts in a low density 
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Fig. 6. B — K vs. if for the if-selected samples. The 
bright if-sele cted s ample (K < 16 mag) is obtained by 
Huang et al. ( 1997 ). The no-evolution and passive evolu- 
tion models are plotted in separate panels. In each panel, 
the upper solid line is the 1.5 L* elliptical galaxy model; 
the middle line is the L» Sbc galaxy model; and the lower 
line is the 0.25 L* Sd/Irr galaxy model. The dashed lines 
link the colors among the different models at z = 0.5 and 
z = 1.0 (for ho = 0.5, qo = 0.05). Though the dashed lines 
are not exactly redshift contours, they roughly outline a 
rcdshift distribution on the color-magnitude relation. The 
cross is a medium B — K color at each magnitude bin. 

universe reasonably well. The no-evolution models are sig- 
nificantly flatter than the observational data. The i?-band 
number counts show the highest excess and may favor an 
even higher normalization. The if-band counts have the 
lowest excess, but the no-evolution model still significantly 
under-predicts the counts for K > 16 mag. 

We note that Metcalfe et al. (1996) and McCracken 
et al. (2000) are able to fit their no-evolution model to 
the observed if-band number counts at faint magnitudes, 
while their passive evolution model predicts more galax- 
ies than observed. This discrepancy between our model 
and theirs is due to the differences in the galaxy spectral 
evolution models in the if-band. 

Although our passively evolving model is able to fit the 
observed number counts in the optical and near-infrared 
bands well, this does not necessarily mean that the pas- 
sive evolution model is a correct physical description for 



Fig. 7. R — K vs. K for the fC-selected samples. The 
models and the medium colors are plotted in the same 
way as in Fig. [(]. 



galaxy evolution at intermediate redshifts. This is rather 
evident in Fig. || where we compare the spectral mix of our 
galaxies, as determined from their SEDs, with the pre- 
dicted type mix from the no-evolution and the passive- 
evolution models. In all three bands this comparison is 
done in the 1.5 mag interval just brighter than the com- 
pleteness limit of our data, where the statistics are best. 
The data and model predictions obviously disagree. The 
passive evolution model predicts that early type galax- 
ies are mainly responsible for the excess over no-evolution 
model at faint magnitudes and high redshifts, Stellar evo- 
lution would make these early type galaxies much brighter 
at high redshifts because their stellar populations were 
younger. However, the spectral type mix of our number 
counts shows that even in the if-band spirals and SD /Irr 
are by far the main contributors to the counts. This can be 
understood if a population of star forming dwarf galaxies 
at z < 1.0 dominates the number counts at faint magni- 
tudes. 

Another way of studying galaxy populations and their 
evolution is through color-magnitude diagrams. Fig. |^ 
(B - K vs. K) and Fig. %{R- K vs. K) are the color 
magnitude diagrams for the if-selected sample. The B — K 
colors for the fC-selected bright sample obtained by Huang 
( 1997 ) are also plotted in Fig. ||. To compare with the data, 
we derived the B — K and R — K colors as a function of 
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Table 6. Colors for E/SO, Sbc and Sd/Irr at z=0 



E/SO Sbc Sd/Irr 
K 4.18 3.86 3.28 
K 2.62 2.63 2.30 



B 
R 



apparent magnitude K for a 1.5 L* elliptical galaxy, a L* 
Sbc galaxy and a 0.25 L* Sd/Irr galaxy and plot them 
together with the observed color-magnitude data. Those 
three galaxy types have been used to calculate the model 
counts presented above, and they represent the popula- 
tions with the reddest, medium, and bluest colors (the 
B — K and R — K colors at z = are given in Table ^J). 
Provided the models accurately represent the color evo- 
lution of galaxies, the E and Sd/Irr model colors should 
bracket the observed distribution of colors and Sbc model 
colors trace the median of the color distribution. As seen in 
Fig. ||, the no-evolution and passive evolution models both 
bracket the B — K color distribution at bright magnitudes 
very well, and both predict a break in the B — K color at 
K ~ 18 mag. For the bright sample of Huang (1997), there 



is no significant difference between the no-evolution and 
passive evolution models. However, at K > 18 mag the 
non evolving Sd/Irr galaxy model colors are redder than 
a significant fraction of the observed colors. The passive 
evolution model predicts much bluer colors for these three 
types of galaxies at K > 18 mag and gives a better repre- 
sentation of the observed color distribution. 

The dashed lines in Figs. || and [t| link the different 
models at z = 0.5 and z — 1.0 (for ho = 0.5, qo = 0.05). 
We can see why even the passive evolution model was un- 
able to predict the type mix in Fig. |[ the ellipticals were 
brighter in the past and thus visible out to higher red- 
shifts [z > 1), making a significant contribution to the 
faint number counts. There are two obvious ways to ex- 
plain this discrepancy between the observations and the 
models: 

(i) Massive galaxies with high star formation rates at z > 1 
appear bluer and are accordingly classified as SED type 
later than E/SO, even though they will end up as early 
type galaxies at z = 0. 

(ii) less massive galaxies do not undergo significant star 
formation until relatively late times [z < 1). 

The good statistics from large galaxy samples and 
the spectral type studies favor the second alternative. 
However, it will be necessary to compare the model spec- 
tral type with the actual morphological Hubble type of the 
K > 18.0 mag galaxies before this question can be settled. 

5. Summary 

We derived the B-, R- and if-band number counts in 
the Calar Alto Deep Imaging Survey. The total coverage 
for the -K'-band imaging is about 0.2 deg 2 . The B-band 
number counts cover the range from B = 16.75 mag to 
B = 24.75 mag; the i?-band number counts cover from 
R = 15.75 mag to R = 23.25 mag; and the if -band num- 



ber counts cover from K — 13.75 mag to K = 19.75 mag. 
Our number counts at optical wavelength match the num- 
ber counts of other surveys very well. Our if -band counts 
differ from the results of previous surveys by up to 50%. 
Because of the large area covered by our survey down to 
K = 19.75 mag, the CADIS number counts in K set a 
new standard in the range 17.75 < K < 19.75 mag, where 
our survey represents the largest if-band survey to date. 
Between K < 16.5 mag and K > 17.5 mag our counts dis- 
play a break in the logarithmic slope from 0.64 to 0.36, 
we estimate this break occurs around K = 16.7 mag. 

In order to get a first insight in the underlying evolu- 
tion process we also adopt no-evolution and passive evo- 
lution models to fit our optical and near-infrared number 
counts. We find that both our optical and near-infrared 
number counts have an excess over the no-evolution model 
and can be simultaneously well fitted by the passive evo- 
lution model. The passive evolution model, however, pre- 
dicts that the excess of the number counts at faint magni- 
tudes is mainly due to the early type galaxies which were 
brighter at high redshifts. This is in contradiction to the 
observed spectral type mix contributions to our counts, 
which seems to be dominated by spirals and star-forming 
galaxies down to the faintest limits of our survey. 
Thus we conclude that both non-evolving and passive evo- 
lution are both in contradiction with our observed number 
counts and the galaxy type mix at the faint magnitude end 
of our survey. The most sensible explanation for this find- 
ing would be a population at z < 1.0 with violent star 
burst events. These make them very prominent at inter- 
mediate redshifts while locally (at z ~ 0) they appear at 
such faint ends of both, the total luminosity and surface 
brightness range that they are underrepresented in current 
redshift surveys. 



References 

Abraham R.G., Tanvir N.R., Santiago B.X., et al, 1996, 

MNRAS, 279, L47 
Arnouts S., de Lapparent V., Mathez G., et al., 1997, A&A 

Supl., 124, 163 

Arnouts S., D'Odorico S., Cristiani S., et al., 1999, A&A, 341, 
641 

Bershady M.A., Lowenthal J.D., Koo D.C., 1998, ApJ, 505, 50 

Bertin E., Arnouts S., 1996, A&A, 117, 393 

Bertelli G., Bressan A., Chiosi C, Fagotto F., Nasi E., 1994, 

A&AS, 106, 275 
Bizenberger P., McCaughrean M., Thompson D., Birk C, 1998, 

SPIE Conference Proceedings, vol. 3354, 109 
Bruzual G., Chariot S. 1993, ApJ, 405, 538 
Casali M.M., Hawarden T., 1992, JCMT UKIRT Newsletter, 

4, 33 

Chariot S., 1996, in From Stars to Galaxies: The Impact of 
Stellar Physics on Galaxy Evolution, ed. C. Leitherer, 
U. Fritze-von-Alvensleben, and J. Huchra, ASP Conference 
Series Vol. 98 

Cowie L.L., Gardner J. P., Hu E.M., et al., 1994, ApJ, 434, 114 
Cowie L.L., Songaila A., Hu E.M. Cohen J.G., 1996, AJ, 112, 
839 



10 



J.-S. Huang et al.: The CADIS near-infrared and optical number counts 



Cowie L.L., Hu E.M., Songaila A., Egami, E., 1997, ApJ, 481, 
L9 

Dickinson M., 2000, in Building Galaxies: From the Primordial 
Universe to the Present, ed. F. Hammer, T.X. Thuan, V. 
Cayatte, B. Guiderdoni, & J. Tranh Than Van, (Paris:Ed. 
Frontieres) 

Djorgovski S., Soifer B.T., Pahre M.A., et al., 1995, ApJ, 438, 
L13 

Ellis R.S., Colless M., Broadhurst T.J., Heyl J.S., Glazebrook 

K., 1996, MNRAS, 280, 235 
Ellis R.S., 1997, ARAA, 35, 389 

Fioc M., Rocca-Volmerange B., 1997, A&A, 326, 950 

Fried J.W., von Kuhlmann B., Meisenheimer W., et al., 2001, 

A&A, accepted (astrph/0012343) 
Gardner J.P., Sharpies R.M., Carrasco B.E., Frenk C.S., 1996, 

MNRAS, 282, LI 
Gardner J.P., 1998. PASP 110,291 

Gardner J. P., Cowie L.L., Wainscoat R.J., 1993, ApJ, 415 L9 
Glazebrook K., Ellis R.S., Santiago B.X., Griffiths R., 1995, 

MNRAS, 275, L19 
Hogg D.W., Pahre M.A., McCarthy J.K., et al., 1997, MNRAS 

288, 404 

Huang J.-S., Cowie L.L., Gardner J.P., et al., 1997, ApJ, 476, 
12 

Huang J.-S., 1997, Ph.D. Thesis, University of Hawaii 
Huang J.-S., Cowie L.L., Luppino G.A., 1998, ApJ, 496, 31 
Kinney A.L., Calzetti D.A., Bohlin R.C., et al., 1996, ApJ, 467, 
38 

Kiimmel M.W., Wagner S.J., 2000, A&A, 353, 867 

Kiimmel M.W., Wagner S.J., 2001, A&A, submitted 

Lilly S.J., Tresse L., Hammer F., Crampton D., LeFevre O., 

1995, ApJ, 455, 108 

McCracken H.J., Metcalfe N., Shanks T., et al., 2000, MNRAS, 
311, 707 

Meisenheimer K., Roser H.-J., 1986, in The Optimization of the 
Use of CCD Detektors in Astronomy, ed. J. -P. Baluteau 
& S. D'Odorico (Munich:ESO-OHP), 227 

Meisenheimer K., Beckwith S., Fockenbrock R., et al., 

1998, in ASP Conf. Ser. 146, in The Young Universe, 
ed S. D'Odorico, A. Fontana, & E. Giallongo (San 
Francisco ASP), 134 

Metcalfe N., Shanks I., Fong R., Roche N., 1995, MNRAS, 273, 
257 

Metcalfe N., Shanks I., Campos R., Fong R., Gardner J. P., 

1996, Nature, 383, 236 

Minezaki T., Kobayashii Y., Yoshii Y., Peterson B.A., 1998, 
ApJ, 494, 111 

Moustakas LA., Davis M., Graham J.R., et al., 1997, ApJ, 
475, 445 

Prandoni I., Wichmann R., da Costa L., et al., 1999, A&A, 
345, 448 

Saracco P., Iovino A., Garilli B., Maccagni D., Chincarini, C, 

1997, AJ, 114, 887 

Saracco P., D'Odorico S., Moorwood A., et al, 1999, A&A, 349, 
751 

Smail I., Hogg D.W., Yan L., Cohen J.G., 1995, ApJ 449, L105 
Songaila A., Cowie L.L., Hu E.M., Gardner J.P., 1994, ApJS, 
94, 461 

Steidel, C. C, Giavalisco, M., Petini, M., Dickinson, M., 

Adelberger, K. L., 1996 ApJ, 462, L17 
Steidel, C. C, Adelberger, K. L., Giavalisco, M., Dickinson, 

M., Petini, M., 1999, ApJ, 519, 1 
Szokoly G.P., Subbarao M.U., Connolly A.J., Mobasher B., 

1999, ApJ, 492, 452 



Thommes E., Meisenheimer, K., Fockenbrock R., et al., 1998, 

MNRAS, 293, L6 
Thompson D., Beckwith S.V.W., Fockenbrock R., et al., 1999, 

ApJ, 523, 100 

van den Bergh S., Abraham R.G., Ellis R.S., et al., 1996, AJ, 
112, 359 

Wainscoat R.J., Cowie L.L., 1992, AJ, 103, 332 
Williams R.E., Blacker B., Dickinson M., et al., 1996, AJ 112, 
1335 

Worthy G., 1994, ApJS, 101, 181 

Wolf C, 1998, Ph.D. Thesis, University of Heidelberg 

Wolf C, Mundt R., Thompson D., et al., 1999, A&A, 343, 399 



